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ABSTRACT 

The 2-ni Liverpool Telescope (LT), owned by Liverpool John Moores University, is located in La 
Palma (Canary Islands) and operates in fully robotic mode. In 2005, the LT began conducting an 
automatic GRB follow-up program. On receiving an automatic GRB alert from a Gamma-Ray Ob- 
servatory {Swift, INTEGRAL, HETE-II, IPN) the LT initiates a special override mode that conducts 
follow-up observations within 2-3 min of the GRB onset. This follow-up procedure begins with an 
initial sequence of short (10-s) exposures acquired through an r' band filter. These images are re- 
duced, analyzed and interpreted automatically using pipeline software developed by our team called 
"LT-TRAP^ (Liverpool Telescope Transient Rapid Analysis Pipeline); the automatic detection and 
successful identification of an unknown and potentially fading optical transient triggers a subsequent 
multi-color imaging sequence. In the case of a candidate brighter than r' = 15, either a polarimetric 
(from 2006) or a spectroscopic observation (from 2007) will be triggered on the LT. If no candidate is 
identified, the telescope continues to obtain z' , r' and i' band imaging with increasingly longer expo- 
sure times. Here we present a detailed description of the LT-TRAP and briefly discuss the illustrative 
case of the afterglow of GRB 050502a, whose automatic identification by the LT just 3 min after the 
GRB, led to the acquisition of the first early-time (<1 hr) multi-color light curve of a GRB afterglow. 

Subject headings: gamma rays: bursts - telescopes - techniques: image processing 
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1. INTRODUCTION 

Since the discovery o f the first optical afterglow 
ijVan Paradiis et al.ll997j) associated with a Gamma-Ray 
Burst (GRB), several robotic telescopes have been de- 
veloped to perform rapid follow-up observations of GRB 
afterglows. Fast robotic follow-ups are now possible be- 
cause of the improved localization of GRBs provided 
by current gamma-ray satellites (HETE-II, INTEGRAL, 
Swift, Interplanetary Network); real-time GRB alerts are 
received from these satel lites through the GRB Coor- 
dinates Network (GCN) (' Barthelmv et alJl2005D within 
a few seconds of the GRB onset and provide positions 
with typical uncertainties of a few arcmin. Observations 
within minutes of the GRB can be crucial to determine 
the properties of the close circumburst environment, po- 
tentially providing clues to the nature of the progenitor. 

To date, only a few GRBs have measurements of their 
optical afterglow starting as early as a few minutes af- 
ter the GRB: 990123, 021004, 021211, 030418, 041006, 
041218, 041219a, 050319, 050401, 050502a, 050525, 
050712, 050713A, 050730, 050801, 050802, 050820A, 
050824, 050904, 050908, 050922C, 051021A (October 
2005). Figure n shows the early afterglow light curves 
for all of them within the first day after the burst. For 
each GRB we collected unfiltered, R, Rc and V observa- 
tions. 

2. THE LIVERPOOL TELESCOPE 
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The Liverpool Telescope (LT), owned and operated by 
Liverpool John Moores University (LJMU), has a 2-m 
diameter primary mirror, altitude-azimuth design, final 
focal ratio f/10. It was designed and built by Telescope 
Technologies Limited, with the robotic control system 
and instrumentation being pr ovided by the LJM U As- 
trophysics Research Institute l|Steele et al. ll200^ . It is 
situated at the Observatorio del Roque de los Mucha- 
chos, on the Canary Island La Palma. Due to the fully 
opening enclosure and to the fast slew rate of 2°/s, the 
telescope can start observations within 1-3 min of the 
receipt of the GCN alert notice. LT has five instrument 
ports: four folded and one straight-through, selected by 
a deployable rotating mirror in the AG Box within 30 s. 
The telescope is equipped with the RAT Cam Optical 
CCD Camera, and with the SupIRGam 1-2.5 /im Cam- 
era. A Prototype Spectrograph and an optical polarime- 
ter (RINGO) based on a design bv lClarke fc Neumaverl 
(2002) are expected to be operating later in 2006. A 
higher resolution, higher sensitivity and generally more 
capable spectrograph (FRODOSpec) is being developed 
for deployment in 2007. Table Q] reports the main char- 
acteristics of each instrument. In addition to the GRB 
follow-up progra m, which has high prior ity on the LT, the 
robotic control (Fraser & Steele '200 3) and a utomated 
scheduler (Ste ele & Carter 1997: Fraser fc Steele .200|) 
allow the LT to conduct an optimized science program of 
observations of time variable sources over a wide range 
of timescales from Target of Opportunity (e.g. novae, 
'WfW)'ef%M?Sfl;^iiHfsWiH,l^es in gravitational lenses) to longer- 
term monitoring (e.g. AGN, variable stars), when no 
GRB observations are required. 
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Guidorzi et al. 
Table 1. Liverpool Telescope Instrumentation 



Instrument 



Description 



RATCam Optical CCD Camera - 



2048x2048 pixels, 0.135" /pixel, 
8 filter selections (u', g', r', i', z 
from LT first light, July 2003 



FOV 4.6'x4.6', 
, B, V, Ha, ND2.0) 



256x256 pixels, 0.4"/pixel, FOV 1.7'xl.7', 
Z, J, H, K' filte rs - from late 2005 
1.7" fibres. 



SupIRCam 1 - 2.5 micron Camera ■ 
(with Imperial College) 
Prototype Spectrograph - 
(with University of Manchester) 



49, 1.7" fibres, 512x512 pixels, R- 
3500 < A < 7000 A- from 2006 



=1000; 



Ring Polarimeter based on the design of ^la^^^^Neuma^e^^OO^ 

- from 2006 



RINGO optical polarimeter 



FRUDUSpec Integral field 
double beam spectrograph - 
(with University of Southampton) 



=4000, 8000; 
4000 < A < 9500 A 



from 2007 




Time since GRB (days) 

Fig. 1. — Early light curves (unfiltered, R and V) for a 
set of GRBs with detections within minutes of the gamma- 
ray event. Big empty upside down triangles show the case of 
GRB 050502a (filter r') robotically detected and followed- up by 
LT (from Guidorzi et al. 1 2005a)). All the measurements shown 
have been taken from GCN circulars (Nos. 248, 257, 307, 715, 
1572, 1589, 1596, 1732, 1738, 1739, 1748, 1752, 1759, 1821, 2773, 
2776, 2784, 2798, 2799, 2821, 2875, 2897, 3118, 3165, 3480, 3485, 
3487, 3489, 3492, 3575, 3584, 3588, 3596, 3601, 3604, 3646, 3704, 
3705, 3711, 3712, 3716, 3717, 3718, 3720, 3723, 3726, 3728, 3733, 
3736, 3739, 3741, 3744, 3745, 3746, 3756, 3765, 3775, 3778, 3829, 
3834, 3836, 3838, 3853, 3863, 3864, 3865, 3868, 3869, 3870, 3880, 
3883, 3896, 3907, 3917, 3929, 3943, 3944, 3945, 3947, 3950, 3953, 
3960, 4011, 4012, 4015, 4016, 4018, 4023, 4026, 4027, 4032, 4041, 
4044, 4046, 4 048, 4049, 4095, 4120, 41 21), except for G RB 030418 
JRvkoff et aU .2004). GRB 041219a IVestrand et al.l .2005.) and 
GRB 05n502a Tnuidorzi et ani2005al) . 



3. REAL-TIME GRB FOLLOW-UP STRATEGY - OVERVIEW 

When a GRB fulfills all the requirements to be ob- 
served with the LT, the GRB over-ride mode is triggered. 
The GRB program takes control of the telescope and 
triggers the so-called "Detection Mode" (SecO, during 
which a small number of Sloan r'-band short-exposure 
images are taken^. These images are then transferred to 
the LT proxy machine, a network-attached on-site com- 
puter. The automatic LT GRB pipeline, called Liver- 
pool Telescope Transient Rapid Analysis Pipeline (here- 
after LT-TRAP), is then invoked to process, analyze and 
interpret these images and returns a value to indicate 
whether any variable optical transient (OT) candidate 

^ From end 2006, this detection mode will be preceded by a short 
30-s exposure taken with the RINGO polarimeter. 



has been identified. The result of this analysis drives the 
automatic selection of the subsequent observing mode as 
follows: 

• identification of a fading OT with r' > 15 ^ Multi- 
Color Imaging Mode fSec. Ifi.ljl : 

• identification of a fading OT with r' < 15 —> Po- 
larimetry Mode (2006)/Spectroscopy Mode (2007) 
(Sec. El; 

• no identification of an OT RIZ Mode fSec. lOji . 

Figure 121 shows a diagram of the overall observing strat- 
egy when the GRB over-ride mode is triggered. 

4. LT GRB OVER-RIDE MODE 

The GRB Target of Opportunity Control Agent (GRB 
TOGA) is a shell script that is invoked in response to a 
GRB alert being received at LT via a socket connection. 
This script uses the TOGA of the LT Robotic Gontrol 
System (RGS) to take control of the telescope, provided 
that the GRB position is observable and other require- 
ments are fulfilled, e.g.: 

• the positional error circle is sufficiently small (error 
radius < 10 arcmin); 

• altitude above horizon is > 20°; 

• the position does not match any known X-ray cat- 
alog source; 

• in the case of Swift alert, the image significance 
should not be low (to-date, no real GRB was found 
to have the low image significance flag raised); 

• none of the five known Soft-Gamma ray Repeaters 
(SGRs), SGR 1900+14, SGR 1806-20, SGR 1627- 
41, SGR 0526-66 and SGR 1801-23, lies in the error 
circle. 

As noted above, the average delay between the arrival 
time of the GCN alert and the start time of observations 
ranges between 1 and 3 min. 

In the case of an SGR lying in the error circle, the GRB 
TOGA triggers a SGR follow-up IR observation with the 
SupIRCam according to the strategy defined in an SGR 
follow-up approved proposal led by the LT GRB group. 
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Currently, the automatic SGR follow-up program is still 
under test, so the presence of an SGR in the error circle 
merely invalidates the GRB follow-up observation. 

When all of the above conditions are fulfilled, the GRB 
TOGA repoints the telescope, then triggers the acquisi- 
tion of the first set of images scheduled in the detection 
mode (following commissioning of RINGO (2006), the 
detection mode exposures will be preceded by an initial 
short polarimetric exposure). The LT-TRAP is then in- 
voked to analyze these images and consequently choose 
the most suitable observing strategy (see Sec. El for a 
detailed description). The RINGO exposure will be an- 
alyzed offline at a later time and plays no part in real- 
time decision making; the motivations of this choice are 
explained in Section 16.21 

The GRB TOGA and the LT-TRAP communicate 
through the noticeboard, which is a plain text file con- 
tinuously updated until the end of observations. This 
not only reports in real time a detailed chronological se- 
quence of all the actions and results acquired through, 
following technical conventions, but also represents a 
useful real-time human-readable interface for the users 
responsible for monitoring the on-going observation se- 
quence. 

5. DETECTION MODE (DM) 

The DM consists of a sequence of two sets of three 10-s 
exposure images taken in r' . These are the first images of 
the GRB field acquired by the telescope. Once each set 
has been acquired and copied to the proxy machine, the 
GRB TOGA invokes the GRB astromctric fit script on 
each image. After object extraction and comparison with 
standard catalogs, an OT candidates list is extracted for 
each set of three images. These OT candidates are then 
processed by the GRB variability script, which aims to 
identify sources common to each image set and evaluate 
possible variability. Finally the DM reports the best OT 
candidate (if any) out of the two merged lists. 

5.1. The GRB Pipeline Astrometric Fit Script 

This shell is responsible for fitting the astrometry and 
identifying possible OT candidates. Below we report the 
detailed sequence of operations performed on each single 
image by this script: 

1. Source extractioi i usin^^ a cu stomized version of 
SExtractor 2.3.2 l)Bert in"199(y) is performed; if the 
number of sources extracted is lower than 3, it does 
not attempt to fit the astrometry. In contrast, if 
more than 20 sources are extracted, only the 20 
brightest ones are used in the astrometric fit pro- 
cedure. 

2. A 9' X 9' portion of the USNO-Bl.O astrometric cat- 
alog (Monet et al. 2003), centered on the centroid 
position of the GRB, is retrieved from a complete 
copy stored locally on the LT proxy machine. In 
case of problems accessing the local catalog, the 
script attempts to retrieve it through the on-line 
interface of the USNO-Bl.O^. 

3. The retrieved catalog image is then cropped to 
match exactly the field of view of the LT. 

http:/ /www. nofs.navy.mil/data/fchpix/ 



4. Optical extinctio n for the field is deter mined from 
the dust maps bv lSchlegel ^TaTI ^f^998^ ■ 

5. An astrometric fit is pe rformed using an adapted 
version of the WGS (Grciscn fc Galabrettal 12002(1 
tool imwcs 3.5.3 LMink ..2002) . The fit is consid- 
ered to have failed if either the number of matched 
sources is lower than 3 or if the fraction of matched 
sources is lower than 50%. 

6. If the astrometric fit is successful, the pixel-to-sky 
coordinates transformation is performe d through 
the WGS tool xy2sky .9. ,5. 7 (|Mink l l2f)0l . based on 
the result of the astrometric fit. A typical value for 
the fit residual is ~0.3 arcsec. 

7. In case of astrometric fit failure, all the image 
sources are considered in principle as potential OT 
candidates and they are reported in the OT candi- 
dates files, although with a lower confidence level 
(see below). The pixel-to-sky coordinate transfor- 
mation is performed assuming the nominal point- 
ing (typical accuracy currently of ^^10-20 arcsec, 
with a target on completion of commissioning of 
the telescope of 2 arcsec). In this case, the script 
ends at this point. 

8. Selection of the OT candidates: if the astrometric 
fit is successful, all the image sources are cross- 
checked with the catalog stars and classified ac- 
cording to a set of criteria. A 17-bit variable is 
assigned to each source and describes its proper- 
ties. Each image source not associated with any 
USNO-Bl.O star is assigned 50 points. Each bit of 
the 17- bit variable is assigned a given number of 
points to be added to the initial 50 points credit 
(see Table 13 . 

9. Selected classes of non-associated image sources 
are a priori excluded from the candidate search. 
Among them, the objects with clearly non-stellar 
PSFs are probably the most important (e.g. cos- 
mic ray events, reflections, failed deblending). A 
mean stellar PSF is defined on each image based 
on real-time associated objects median statistics. 

10. Each selected OT candidate is assigned a confi- 
dence level directly related to the total number of 
points accrued. 

11. Each selected OT candidate is cross-checked with 
the on-line GSG2.3 and 2MASS catalogs and its 
confidence level is adjusted accordingly (see Ta- 
ble 0). 

When the astrometric fit is successful, magnitudes are 
estimated from the USNO-Bl.O stars that have been 
matched. We use as zero point the median value of the 
offset needed to have a match between the instrumental 
magnitude and the R magnitude reported in the USNO- 
Bl.O catalog. The error associated with the zero-point 
is then evaluated in terms of the clipped (2ct) root mean 
square deviation with respect to the USNO-Bl.O. 
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Table 2. LT-TRAP: 17-bit variable assigned to each OT candidate 



Bit n 


2" 


Description 


Script 


Points 


\ 


\ 


LliC CA.Lrl dt, Levi llldglll L Lltlt; lo UidOtJLI. IJy LJl IgliL llci^ilDiJ Lli o 


iT.i5 LI Ullie LI IC 


Q 


2 


2 


ttiG ob j Gct h&s bGGn dGblciidGci 


Astrometric 


-10 


3 


4 




JIO Li UlllC LI IL, 





4 


8 


object trunccited (close to the inicige border^ 


Astrometric 


-50 


5 


16 


ctperture dS/tSy incomplete or corrupted 


Astrometric 


-10 


6 


32 


isoptiota,l da,ta, incompletG or corrupted 


Astrometric 


-10 


7 


64 


memory overflow occurred during deblending 


Astrometric 


-10 


8 


128 


memory overflow during extraction 


Astrometric 


-20 


9 


256 


wasn't possible to check GSC and/or 2MASS 


Astrometric 


-20 


10 


512 


a GSC2.3 object lies close (< 2 arcsec) 


Astrometric 


-50 


11 


1024 


a 2MASS object lies close (< 2 arcsec) 


Astrometric 


-30 


12 


2048 


a bright USNO star {R<13) is close (<50px) 


Astrometric 


-20 


13 


4096 


ASTROMETRIC FIT FAILED Assuming nominal pointing 


Astrometric 


-50 


14 


8192 


Object detected in at least 2 OT lists 


Variability 


-f30 


15 


16384 


Object NOT detected in at least 1 OT list 


Variability 


-10 


16 


32768 


Variability detected over 2—3 images 


Variability 


-1-50 


17 


65536 


Astrometric fit warning 


Variability 






5.2. The GRB Pipeline Variability Script 

After the LT-TRAP script described above has run on 
the three images, the three resulting OT candidates' files, 
each reporting a list of potential OTs found in each sin- 
gle image, are cross-checked through the so-called "vari- 
ability script" . This shell script is invoked by the GRB 
TOGA and aims to recognize objects common to each of 
the OT candidates files, as well as to evaluate possible 
significant magnitude variations between the three im- 
ages. This variation is evaluated by performing a test 
on the zero hypothesis that there is no significant vari- 
ation. If the three values of magnitude are found to be 
consistent with a constant value with a significance lower 
than 1%, the variability bit is raised (Tabled). Each OT 
candidate confidence level is finally adjusted on the basis 
of the variability results (bit number 16, see Table [5J). 

Table 12 reports the meaning of each bit of the 17-bit 
variable assigned to each OT candidate and the points 
assigned to each flag. The confidence level is evaluated 
by dividing the total points by 100: when this ratio is 
< 0, the confidence level is set to zero; when the ratio 
is > 1, the level is assigned 1.0 (highest confidence level 
possible). The final threshold on the confidence level 
for an OT candidate to trigger the Multi-Color Imaging 
Mode (Sec. 16.1(1 is currently set at 0.7. In case of more 
than one OT candidate with confidence level above the 
threshold, the GRB TOGA accepts that with the highest 
value. If the ambiguity is still unresolved, it assumes the 
first on the list. However, all the information regarding 
the other candidates is reported in the noticeboard file 
as well as in the product files. 

The meaning of each single bit described in Table El is 
self-explanatory, apart from the last bit, called "Astro- 
metric fit warning" , which requires an explanation. This 
bit is raised whenever at least one of the three images 
has a failure in the astrometric fit and at least another 
one has a successful astrometric fit. In this case, the 
identification of the same object in different images is 
performed on the basis of pixel coordinates, while the 
sky coordinates are taken from the image with the cor- 
rect astrometry. However, so far this case has turned out 
to be very rare: either all or none of them are success- 



fully astrometrically fitted. For this reason, the current 
points assigned to this bit are zero. 

To test the sensitivity of the script, we simulated on 
real images an OT fading according to a power law, 
F (x {t — to)~" (where to is the GRB onset time). The 
script detected significant variability for a ~l-2, for 
objects at least as bright as r' ~16-17, provided that 
the DM is triggered at [t - to) < 200 s. Also, in the 
case of GRB 050502a, where the DM was triggered at 
t—to = 187 s, the LT'-T'i?v4P detected significant variabil- 
ity: the magnitudes evaluated by the LT-TRAP were for 
the three DM images: r[ = 15.57±0.03,r^ = 15.66±0.03 
and Tg = 15.83 ± 0.03. The corresponding test gave 
the following result: x^/dof = 36.2/2, with a significance 
of 1.4 X 10"*^. 

It must be pointed out that the systematic error af- 
fecting the USNO-Bl.O magnitudes is around 0.3 mag; 
furthermore, it is not possible to correct for color terms at 
this stage. However, this is not relevant, as the purpose 
of the script is to establish relative variations, regardless 
of the absolute values. To account for possible varia- 
tions in the zero point, the script adopts the following 
approach: it calculates the root mean square deviation 
of the three zero points and adds it in quadrature to the 
statistical error of each of the three magnitudes. 

6. OBSERVATION MODES 

After the DM has completed its imaging and analysis 
sequence, the appropriate observation mode is triggered 
automatically on the basis of the results of the analy- 
sis performed by the LT-TRAP on the DM images (see 
Fig. El and Table EJ. 

6.1. Multi-Color Imaging Mode (MCIM) 

This mode is triggered by the GRB TOGA when the 
DM finds an OT candidate with the required confidence 
level. It consists of a multi-color imaging sequence cy- 
cling through three Sloan g' , i' and r' filters until the 
end of the first hour. The ordered sequence, g', i' and 
r', is repeated with different exposure times: the first 
two times with 30-s exposure, then it continues with 60- 
s exposure cycles up to 18 times. Before May 16, 2005, 
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Telescope over-ride 













^ Short RINGO Exposure (2006) 



RIZ Mode 
(RIZM) 



LT 

FTN 

LT & FTN 




Multi-Color 
Imaging Mode 

(MCIM) 



Polarimetry Mode (PM) (2006) 
Spectroscopy Mode (SM) (2007) 



z',r',i' (B,R,i') Light Curves 
or deep Upper Limits 



1 



Spectroscopy 
Mode (SM) 



Polarimetry 



1 



Spectra 



Multi-Color 
Light Curves 



Fig. 2. — Diagram of the LT- T_Ry4P observing strategy for both LT and FTN facilities. 



this mode operated in the same way but used a slightly 
different se t of fihcrs, BVr'i', as used for GRB 0505G2a 
ijGuidorz i gt al. 2()Q^^ ) . The motivation for changing and 
optimizing the filter choice (post May 16) was threefold: 
1) the greater sensitivity at g' band, 2) the small separa- 
tion between B and V filter central wavelengths (444.8 A 
and 550.5 A, respectively) and 3) the recognized need for 
denser sampling of the light curves at early times. 

6.2. Polarimetry Mode (PM) 

A key unsolved problem in the astrophysics of GRBs 
is the geometry of the ejecta and the role of magnetic 
fields. Direct measurements of polarization at early times 
provide stringent constraints on current jet and internal 
shock models. As described in Section 12 we will have 
unique access to a simple optical polarimeter, RINGO, 
on the LT from the end of 2006. Polarimetry Mode 
(PM) will be triggered following the identification by 
LT-TRAP of a candidate OT with r' < 15 mag and 
high significance, and will consist of an automatic series 
of polarimetric exposures taken with RINGO, with in- 
creasing exposure times: 30s, 60s, 120s, 240s, etc. These 
exposures will not be reduced and analyzed in real-time, 
but post-processing will extract the r' band light curve 
and any polarization evolution of the OT. Although less 
sensitive than RATGam, RINGO will detect polarization 
to a few % in r' < 15 mag OTs in a 30-s exposure. In 
combination with the initial polarimetric exposure that 
will be obtained immediately before the DM begins, this 



mode aims to measure the polarization properties from 
the first minutes after the burst. 

6.3. Spectroscopy Mode (SM) 

Following commissioning of FRODOSpec (2007), the 
identification of a candidate OT with r' < 15 mag and 
high significance will result in the triggering of an au- 
tomatic spectroscopic observation. The goal of the ini- 
tial spectroscopic observations will be to obtain an early 
redshift estimate, whilst subsequent spectra taken dur- 
ing the first hour will provide valuable information on 
the time evolution of the GRB environment and light 
curve. Short exposure r' will be interspersed with spec- 
troscopic observations to ensure the OT has not faded 
below the sensitivity limit of the spectrograph. In the 
case of a rapidly fading OT, the observing mode will 
revert to MCIM after the initial spectrum has been ob- 
tained. 

6.4. RIZ Mode (RIZM) 

When the DM does not identify an OT candidate, the 
GRB TOGA invokes the "RIZ Mode" . This consists of a 
sequence of 2-min exposure r', i' and z' images: first, one 
r', one i' and one z', then two r' followed by two i' and by 
two z', then three r' followed by three i' and three z' and 
so on. Prior to November 2005, this mode consisted of 
r'i' cycles only; we have subsequently adapted this mode 
to include additional z' exposures that are inserted at the 
end of each r'i' cycle to enable identification of possible 
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Table 3. LT-TRAP Strategy: observing modes for different facilities 



Observing Liverpool Telescope Faulkes Telescope North 

Mode 



DM 2 X (3x10 s in r') (3x10 s in R) + (1x10 s in B, V, i') 

RIZM^ (n X 2 min) in (r', 2')=^, n = 1, 2,3,. . . 

MCIM 2 X (30 s in g',r',i') + n x (1 min in g',r',i') cycles of B,R,i': lx30s, Ixlmin, nx (2 min, 3 min) 

SM For OTs with r' < 15 mag (expected in 2007) R < 15 mag 

PM 30s, 60s, 120s... for OTs with r' < 15 mag (expected in 2006) - 



''Original RIM amended to RIZM by inserting additional 2' exposures at the end of each r'i' cycle (from November 2005). 



higher redshift bursts (z < 6) that may be faint in r' but 
brighter in z' . 

This mode is designed to be sensitive to possible faint 
optical afterglows (whether intrinsically faint or at higher 
redshift), allowing r',i',z' band light curves to be con- 
structed when the OT is detectable in longer exposures 
but too faint to have been identified automatically by 
the detection mode, or to provide deep upper limits at 
early time when no OT is present. Under average condi- 
tions, i.e. seeing of 1 arcsec, 2-min exposure limiting 
magnitudes are around r' ~21-22, i' ~20-21 and z' ~20. 

7. LT-TRAP DEPLOYED ON OTHER FACILITIES 

The LT is the prime member of the global 2-m robotic 
telescopes network called "RoboNet-l.ff^ This network 
includes two other facilities: the Faulkes Telescope North 
(FTN) located in Maui (Hawaii) and the Faulkes Tele- 
scope South (FTS) in Siding Spring (Australia), both 
mainly supported by the Dill Faulkes Educational Trust 
and intended for use by UK schools and various scientific 
observing programs. All the three telescopes are oper- 
ated remotely through the Operations Management Cen- 
tre in Liverpool JMU and, when fully commissioned, will 
all operate in fully robotic mode. Under the RohoNet-1.0 
project, comprising a consortium of ten UK universities 
funded by the UK PPARC, a fraction of time on the 
two FTs is reserved for pos t first hour GRB foUqw-up by 
consortium members (see IGomboc et al.l (j2004l I2005a|l 
for details). 

Following the successful operation of LT-TRAP on the 
LT, it was deployed on the Faulkes Telescopes, with a 
slightly-modified observing mode choice following the OT 
identification stage (see Fig. \^ . FTS is currently under- 
going commissioning so is not yet available robotically; 
the FTN is fully robotic and has performed a number of 
follow-up observations minutes after a GRB; for this rea- 
son, hereafter we will limit our discussion to FTN. The 
rapid (< 1 hr) robotic follow-up observations of GRBs on 
FTN is a collaborative project between Liverpool John 
Moores University and the University of Leicester. 

The optical filters available on the FT are different 
from the LT: Bessell BVR, Sloan i' , neutral U, Olll and 
Ha. The DM consists of only one set of three 10-s im- 
ages in R. Regardless of the results of the DM, three 
more 10-s images are acquired in the following filters, re- 
spectively: B, V, i' . In the case of identification of an 
OT candidate with R < 15, the GRB TOGA wiU trig- 
gers a spectroscopic observation. As in the case of LT, 

* See |http://www.astro.livjm.ac.u k/RoboNet| 



the threshold on the confidence level for an OT candi- 
date to be considered as a good candidate is currently 
set to 0.7. The logic of the selection of the OT candidate 
is exactly the same as for the LT (Sec. 15.2(1 . Presently, 
regardless of the possible identification of an OT candi- 
date with i? > 15 and confidence level > 0.7 from the 
DM, the GRB TOGA triggers a sequence of multi-color 
observations according to the following scheme: cycles of 
BRi' with increasing exposure times: 30 s (once), 60 s 
(once), then it goes on continuously switching between 
120 s and 180 s exposures in order to limit the effects of 
tracking problems. 



8. ESTIMATION OF THE LT-TR^P RELIABILITY 

Although the number of GRBs so far observed is still 
small, we have attempted to estimate the reliability of the 
OT candidates automatically identified by the LT-TRAP 
as follows. For each of the GRBs followed up robotically 
within minutes, we counted the OT candidates identified 
as a function of their confidence level. When the confi- 
dence level is above the threshold (currently set to 0.7), 
the OT is taken to be a true OT candidate (TOG), oth- 
erwise it is a fake OT candidate (FOG). When a true af- 
terglow is present, we call it a true OT (TO); similarly, if 
any non-GRB source that could mimic a GRB afterglow 
is present, we call it a fake OT (FO). Let pT be the prob- 
ability that a TO is automatically identified as a TOG. 
The closer to 1, the better the LT-TiJ^P capability of 
identifying genuine GRB afterglows. Likewise, let pp be 
the probability that a FO is automatically identified as 
a FOG. The closer to 1, the better the LT-TRAP ca- 
pability of rejecting fake GRB afterglows. In practice we 
do know that and are somehow correlated: a low 
threshold on the confidence level would increase pt but 
also decrease pp. The opposite is true for a high value of 
the same threshold. So far, the current value (0.7) has 
turned out to be a good trade-off. However, as a first ap- 
proximation we treat them as unrelated and eventually 
we discuss the opposite case of total correlation. Let rix 
and Tot the total number of TO observed and the total 
number of TO identified as TOG, respectively. Likewise, 
let np and mp the total number of FO observed and the 
total number of FO identified as FOG, respectively. The 
probability Pt('^t, "t^tIpt) iPFinp,mp\pp)) to identify 
TOT (mp) TOG (FOG) out of ut TO (np FO), given 
Pt (pf), is given by the binomial distribution (under the 
sensible assumption that different GRBs are independent 
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from each other): 

Pi{ni,mi\pi) = 



We do not know either or pp, but so far we counted 
the following: rtx — 2, = 2; np = 15, mp = 13 (in 
fact LT-TRAP so far correctly identified both r' < 19 
OTs and failed to reject 2 out of 15 fake OTs). From the 
Bayes theorem, we derive the probability Pi{pi\ni, mi) to 
have Pi , given rii and mi : 



Pi{pi\ni,mi) 



Pi{ni,mi\pi) P{pi 
P{ni, mi) 



,(* = T,F) (2) 



We then apply the maximum likelihood method and find 
the value for pi which maximizes Pi{pi\ni,mi). To this 
end, the denominator of eq. Elcan be ignored. Concern- 
ing the knowledge of the prior P{pi), in both cases as 
a first guess we may assume the most simple case of a 
uniform distribution. Then we refine it thanks to the 
information derived from our tests, according to which 
a sensible approximation is given by a truncated nor- 
mal distribution centered on 0.8 and sigma of 0.3. How- 
ever, it must be pointed out that the prior strongly de- 
pends on the kind of afterglow population one considers. 
For this reason, we consider these two opposite cases of 
uniform and truncated normal distributions. When a 
uniform prior is assumed, the best value is straightfor- 
wardly given hy Pi = mi /nil px = 2/2 = 1 and pp — 
13/15 ~ 0.87. When we assume a more refined prior, 
i.e. P{pi) ~ A^(0.8,0.3), the best values are pr ^ 0.98 
(Fig. 131 solid line) and pp — 0.86. Eventually, we can 
provide confidence intervals on both p^ and pf at given 
levels from the distributions obtained with the truncated 
normal priors: pt = OM^oTe and pf = OM^o fe (50% 
CL); PT = 0.981°:!^^ and pp = OMto fr (90% CL). 

If we take into account the correlation between p^ and 
Pf, the total probability of eq.Elis replaced by the more 
general one: 

PT,FiPT,PF\nT,m'T,nF,mF) = 

PT{nT,mT\pT) PF{nF,mF\pF) P{pt,Pf) 



PTinT,mT,nF,mF) 



(3) 



At this stage, it is difficult to model the anticorrelation 
between px and pf in the bivariate prior P(px,Pf), par- 
ticularly for the small sample of GRBs so far collected. 
We then studied the opposite case of a total anticorrela- 
tion between the two variables, in which pF is completely 
determined by the value of px, by solving eq. ^ 



/■I rPT 

/ PF{p'\nF,mF)dp' - / Pt{p'\ 

JpF Jo 



rix, TOx) dp' (4) 



where Pi{p'\ni,mi) is the same expression derived in 
eq. 13 assuming truncated normal priors. The mean- 
ing of cq. 2]is that, for each value p'rp, p'f is determined 
so that the probability that pF > Pp, is equal to the 
probability that px < p't', it is straightforward to ver- 
ify that the greater p'p, the smaller p!^ and vice versa. 
This is just one sensible assumption of total anticorrela- 
tion among other possible alternative assumptions. Fig- 
ure 131 shows the probability density for px derived in 
this case of total anticorrelation (dashed line). The best 




Fig. 3. — Probability density function of px assuming a truncated 
normal prior, A'^(0.8, 0.3), in two cases: px a-nd Pf uncorrelated 
(solid line) and totally correlated (dashed line) (see text). 



values for px and pf and their confidence intervals de- 
rived in this case (truncated normal priors) are the fol- 
lowing: Pt = 0.831° °^ and pp = 0.83t°;°^ (50% CL); 

PT = 0.83l°;^^ and pf = 0.83t°;°^ (90% CL). 

Remarkably, in spite of the opposite cases of variables 
Pt and pp either totally unrelated or totally anticorre- 
lated, the results seem to give consistent confidence in- 
tervals, as shown by Fig. 13 The wide confidence interval 
of Px is understandably due to the low number of TO 
observed so far, rather than the limits of the LT-TRAP. 
We expect to give more precise estimates of both proba- 
bilities as soon as the sample of GRB afterglows detected 
is more numerous. 

9. SOME RESULTS AND FUTURE PLANS 

To-date (October 2005), over a period of two months, 
LT followed up in a fully automated fashion using the 
real-time pipeline six GRBs and FTN followed up five 
more. Table 01 reports for each of these eleven follow- 
up observations the start time from the GRB trigger, 
the filters used, the number of frames acquired and 
the GCN circulars issued by the LT/FTN collaboration 
reporting on the corresponding LT/FTN observations. 
The most successful case so far is represented by the 
automatic detection o f the afterglow of GRB 050502a 
fGomboc et al.ll2005 b') followed by the acquisition of the 
iirst early (< Ihr) multi-color light curve (Guidorzi ct a^ 
[2005a|) . This latter case is also notable for longer term 
follow-up using the RoboNet-1.0 netwo rk. The other case 
of pro mpt detection is GRB 050713a l|Monfardini et alJ 
l2005aj) . although in this case the afterglow was not iden- 
tified automatically because of the poor quality of the 
images affected by the presence of a ^ ~ 6 star in the 
field, but recognized afterward by visual inspection. 
The power and robustness of the LT-TRAP is illus trated 
in the case of GRB 050730 (|Gomboc et al.ll2n05d) : this 
GRB occurred before twilight at the LT site, so was 
triggered manually and, although the OT was not rec- 
ognized during the DM because of the high sky back- 
ground due to the sky not yet being completely dark, it 
was automatically identified in the subsequent RIM as 
an unknown fading source of r' = 17.3 at t ~ 50 min 
after the GRB. In the other cases, the combination of 
sensitivity and rapid response made it possible to derive 
deep upper limits, as in the cases of GRB 050504 with 
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Table 4. Results of the robotic follow-up observations with Liver- 
pool AND Faulkes North Telescopes (October 2005) 



Telescope GRB Spacecraft Filters Start Result Frames LT/FTN GCN 

(mill) Circulars 



FTN 


050412 


Swift 


BVRi' 


2.5 


R > 18.7 


23 


none 


LT 


050502a 


INTEGRAL 


BVr'i' 


3.1 


r' ~ 15.8 


25 


3325 


FTN 


050504 


INTEGRAL 


BVRi' 


3.7 


R> 19 


38 


3351 


LT 


050520 


INTEGRAL 


r'i' 


4.5 


r' > 16.6 


18 


3437 


LT 


050528 


Swift 


r'i' 


2.5 


r' > 17.2 


36 


3497 


LT 


050713a 


Swift 


r' 


2.4 


r' ~ 19.2 


3 


3588 


FTN 


050713b 


Swift 


R 


3.3 


R > 18.2 


6 


3592 


FTN 


050716 


Swift 


BVRi' 


3.8 


R > 19.8 


20 


3625 


LT 


050730 


Swift 


r'i' 


50 


r' ~ 17.3 


36 


3706 


LT 


050904 


Swift 


r'i' 


3.8 


b 


31 


none 


FTN 


050925 


Swift 


BVRi' 


3.3 


R > 19.0 


12 


4035 



''This corresponds to the time delay with respect to the GRB trigger time. 

''The afterglow of GRB050904 IHaislip et ani2005bl) was found to lie outside the field of view of the LT (3.9 arcmin away from the BAT 
position). 



i? > 19 (IMonfardini et al.ll2005bD and GRB 050716 with 
R > 19.8 l|Guidorzi et alJl200'5b|r 3-4 min after the GRB 
trigger time. 

In the future, we plan to refine the strategy when an 
OT candidate with i? > 15 and confidence level > 0.7 is 
found from the DM; in particular, exposure times will be 
tuned dynamically on the basis of the estimated magni- 
tude of the best OT candidate. Automatic polarimetry, 
near-infrared and spectroscopy follow-up modes will be 
implemented on the LT following completion of instru- 
ment commissioning. 
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